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ABSTRACT 


A computer model was developed to simulate a C.R.P. propeller 
System. A set of simultaneous differential equations was solved 
using CSMP III-the IBM simulation language. The computer model was 
verified against manufacturer's data. Power piston displacement, 
power piston velocity and chamber pressure were determined for a 
nominal supply pressure of 600 PSI and various loadings. Overall, 
the system model and data should prove useful for investigations of 


design variations. 
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iP LN PRObUC ELON 


During recent years, naval applications of gas turbine engine 
systems have become more prominent. At first, the gas turbine engine 
was developed as an aircraft power source. Presently, it is entering 
service as a marine power source. Gas turbine engine EyStenS have 
Several advantages compared to other power systems. The main advan- 
tages are its relative compactness and the fact that the rotating out- 
put of the turbine is usually what is wanted to turn a wheel, propeller, 
generator, or pump. Unfortunately, compared to other systems, there 
is a problem of matching the characteristics of the load and the engine. 
In this respect, controllable reversible pitch (C.R.P.) propellers are 
useful. 

Controllable pitch propellers, unlike fixed pitch propellers, have 
the ability to reverse thrust through a reversal of pitch rather than 
a reversal of rotation. The controllable pitch propeller is much more 
complex than a solid propeiler and contains many more parts [1]. How- 
ever, the operation of a C.R.P. system provides a greater degree of 
flexibility in the maneuvering and docking of ships. 

Since there is an increasing interest in the application of con- 
trollable reversible pitch propellers to naval vessels, the objectives 
of this study were as follows: 

1. to determine an overall analytical model of an operational 

Gan. fa system: 


2. to check the system response for design variations. 


1k 





If. BACKGROUND 


Pee PASIC CONTROL CONCEPTS 

In a control system in which the output appears as mechanical energy, 
Pevewguantity to be controlled can be either load, torque, position, or 
velocity and can be controlled under either open- or closed-loop condi- 
tions. In an open-loop system [2], the input is fed to the controller 
which then operates on the output. The input, being an independent vari- 
able, iS incapable of correcting any errors between the desired and 
Beeval Output. 

The accuracy of an open-loop system is, therefore, greatly affected 
by external disturbances and variations of system parameters. In a closed- 
loop system [2], the output is fed back and compared with the desired 
value. The error between the desired and actual output is fed into the 
control elements which in turn operate on the output, tending to reduce 
or null the error. This is the concept of negative feedback. A block 


diagram of a closed-loop system is shown in Figure l. 
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Figure 1. Block Diagram of a Closed-Loop system. 
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In a position control system in which the manipulated variables are 
pressure and flow, the hydraulic quantities are converted into a dis- 
placement by action on the piston rod or output shaft. The two proper- 
ties of flow which can be manipulated are its direction and magnitude. 
Directional control of flow is used for simple position control systems 


under both open=- and closed-loop conditions [3]. 
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Figure 2. Valve Controlled Piston. 


Dixectional control valves that could be used here are of the three- 
or four-way type. A schematic arrangement of a four-way valve control- 
ling a symmetrical piston is shown in Figure 2. 

When the control valve as shown is in the neutral position, the 
piston is hydraulically locked. A displacement of the valve spool in 


either direction causes the piston to move in the opposite sense. 


ies 





An idealized graph of a position control system which is operated 
under open-loop conditions, by a directional control valve is shown in 


Figure 3. 


OUTPUT 


SNELL 





Figure 3. Position Control System Response (Open-Loop). 


When the valve is displaced from its neutral position, the actuator 
moves at constant velocity in the required direction until the valve 
is returned to neutral. An open-loop system of this type is only suit- 
able for simple applications which do not require accurate piston control. 

Closed-loop position control systems employ controllers having non- 
continuous or "on-off" characteristics. When the actuating signal 
reaches a threshold value, the directional control valve becomes opera- 
tive and causes the actuator to move at constant velocity in the chosen 
direction. It continues to move until the error signal falls below the 
threshold value. This causes the valve to return to its closed position 


and stop the motion of the piston. 
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The most widely used type of control for hydraulic positioning servo 
systems is the control flow type, where both flow direction and magnitude 
are manipulated. The two methods of achieving flow control are: 

1. controlling the displacement of a variable delivery pump; 

2eeeenroteling £lurd through the ree orifices of a control 

valve. 

The two basic differences between pump eee and valve operated 
servo systems can be described as follows: 

1. In pump operated servos, working pressure is variable, 

depending on the external loading, while valve controlled 
Servos Operate at constant system pressure. 
2. In valve controlled servos, power sources are outside 


the loop. 


In this study, a valve controlled servo system was studied. 


B. ANALYTICAL APPROACH 

This anaiytical investigation was directed toward predicting various 
C.R.P. system pressures and motions. An analytical model was developed 
uSing first-order differential equations. The mathematical model, which 
was programmed on a digital computer, predicted the position and velocity 
of the power piston and the pressure at various points in the system as 
Emeeumction of ee Further analysis was performed using the model to 
obtain preliminary predictions of system performance under a variety of 
external loadings. 

To understand the derivation in the following sections, Figure 4 
indicates the geometry of a basic controllable reversible pitch (C.R.P.) 


Propeller system. Beginning with the input side, a force balance on the 
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input servovalve yielded an equation for the load pressure on the input 
piston in terms of the valve motion. The input piston Re moved by load 
pressure, Poa: The valve rod control assembly was moved by the force 
exerted on the input piston. This caused the primary supply pressure, Poe 
to pass through the hub valve into the power piston and move the propeller 
crank-slot pin. This power piston control system, as such, is similar in 
operation to a vehicular power-assisted steering mechanism. 

To facilitate the digital simulation it was assumed that the input 
valve spool displacement, a and the arbitrary load, Bos were step 
inputs. In the process of solution, the system state equations were 


solved using CSMP III, the IBM simulation languace. 
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Tif. CONSTITUENT RELATIONSHIPS 


It was a difficult problem to identify the dynamic characteristics 
of the controllable pitch propeller system because of its essentially 
nonlinear behavior. Identification was partially achieved by assuming 
that these responses were linear and then applying parameter estimation 


for linear systems. 


A. INPUT SERVOVALVE AND PISTON 

The combination of a servovalve and a piston is a common hydraulic 
power element [3]. In analyzing this simple system, three basic rela- 
tionships surface. The first concerns the valve, the second concerns 
the piston dynamics and the third concerns the flow of hydraulic fluid 
between the elements. 


Regarding the valve, the basic equation that pertains is 


208 7 Kaa lee Kea Pea ) 





where K = 70 w yvP 
qa a sa 
7 ee 
ca P 
sa 


Regarding the piston, the basic equation governing value dynamics is 


ae eee 2 
. (2) 


is 





The sum of the force acting on the piston is 


= = Fog - oe = Na x - eee x 3) 
where Soft =A Poa 
Poa = Ke (X - Y) = 0.86 ice - Py) && - Y) 


Combining equations (2) and (3) gives 


A P, =M X+B_ X+K =X+P (4) 
p a V pa pa a 


The flow relationship between the two elements is developed by 


applying the continuity equation to each of the piston chambers. In 


' 


final form, this continuity relation becomes 


V 


° a . 
en ea ee 8.” Pua ) 


h c a icie 
where olla 1s an overall leakage coefficient. 


Equations (1), (4) and (5) are the basic equations of this sub- 
system. 
Bs POWER PISTON 

From the characteristic equation of the critical center spool valve 


(hub valve), one can derive 


oo = kK (<= i) = K om P a (6) 


where R= l40%~ vo 
qm m Ss 
K = 2s 
cm Pp 
Ss 


A, 





In a manner similar to the development of equation (5), the continuity 


expressions for the hub valve can be combined to yield 


V 


=) P (7) 


QO =C A 5. a 


“ma tpkm 





P +A Y + ( 
m m 
where © elk is an overall leakage coefficient. 

The final equation arises by applying Newton's second law to the 


power piston [4]. The resulting force equation is 


= + oe + = + 
Rem Mo ag Bom (X Ys) Kom (X XY) Po (8) 
Equations (6), (7) and (8) are the basic equations of the power 


piston system. 


ee OVERALL SYSTEM 

In analyzing the overall system, two modes of operation became 
apparent. First, if the system input was a forward valve displacement, 
Bees the output was a power piston position and velocity. The second 
mode occurred during recovery from an external loading. For this mode, 
the input was an external load, Fo. At this point, it was necessary 
to simplify the governing equations. Usually, the loading on a valve 


controllec piston system is simpler than the most general in two 


respects [3]. First, significant spring loads are usually absent and 


hence pea. = eri = 0. Second, damping loads are likewise negligible and 
then B = B = Q. 
pa pm 
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Thus, equation (8) reduced to 


ee IL 
Y= ym Pw Fy) (9) 


A combination of equations (6) and (7) produced 








4 8 
P =———(-K  P-A Y+K_X-K_ Y) (10) 
m Mie cem m m qm qm 
Equation (4) was simplified to 
eee ny a) p (11) 
= M fe a p Qa? 
Finally equations (1) and (5) yielded 
, 468. 
Pua i vie Ss ea Poa : A a Ge ee oe 


Equations (9-12) were the governing equations for the C.R.P. system. 
The following is a listing of the constants used in the derivation of 
these governing equations and in the accompanying computer solution. 


The valves were derived from manufacturer's data [5]. 
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A, (AM) 


K oq (CEA) 


ee (VA) 


TABLE I 


BLST OF CONSTANTS 


cross-sectional area of input piston 
cross~sectional area of power piston 
flow gain of input servovalve 

flow gain of power piston 

total volume of input piston 

total volume of power piston 

valve gradient of input servovalve 
valve gradient of power piston 
supply pressure (power piston) 
supply pressure (input servovalve) 
bulk modulus 


total mass of power piston and 
propeller 


total mass of input piston and 
valve rod assembly 


total flow pressure coefficient of 
input piston and servovalve 


total flow pressure coefficient of 
power piston 
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B7.67 in 
LC 2 2 0) an 
aS 
1.54 in /sec/PSI 
Palla ane /see/ PSI 
729671 ae 
Pe 
4326.08 1m 
ee 
a 72) ee Wana 
tee eee 
Lie 93S an an 
400.0 PSI 
600.0 PSI 
100000.0 PSI 


69.8757 Meee 
Zoe 
4.9 lb-sec /in 
gute 
2057.08 in /sec/in 


33426.54 in?/sec/in 





IV. SYSTEM STATE EQUATIONS 


The system equations were first written in a state variable format 


[6]. Thus, in matrix notation, the system state equations became 


X(t) = A X(t) + B U(t) 


Here X(t) was the matrix of system state variables and U(t) was the 
matrix of system input. From the overall system equations (9-12), Y(t), 
Poi (t), X(t) and Pt?) were chosen as the basic problem variables, and 
Xs and Fo became the input variables. Thus, the applicable state 


Variables were 


Y = X(1) X = X(4) P = X(3) 7 X (6) 
We= X(2) X = X(5) ES = X(3) Pos = X(6) 
Y = X(2) x = X(5) 
0) = Aya a 
For example from equation (9) one can write 
An alt 
02) eee = Gin) 
Mo Mo 2 


In a similar way, the system equations were rewritten in matrix 
notation (Figure 5). The accompanying signal flow graph for the over- 
all system equations is presented in Appendix A. The state equations 
were solved for various servovalve stroke values, Xa! and for various 
loading, Pos using the IBM simulation language CSMP III as shown in 


the next section. 
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V. ANALYSIS OF SYSTEM PERFORMANCE 


Manufacturer's data was consulted in order to formulate a power 
piston velocity specification. The maximum power piston velocity 
required was a value of 0.5 in/sec. Common operational velocities 


ranged between 0.4 - 0.5 in/sec. 


meee PROCEDURE 
ie iehethats 
Electro hydraulic servovalves have strokes ranging from 0.005 
mero in. £0r small valves up to 10 gpm capacity to 0.015 to 0.030 in. 
for 50 gpm valves. In this case, a 15 gpm valve was used and input, 
Xa! was set from 0.005 to 0.012 in. 
An external loading, Pos was assumed to range from 0.0 to 
40,000. 1lbf in the simulations. This range of loading could be expected 
in operational conditions. 
ee Outputs 
The scheduling of control parameters is shown in Figure 6. 
Y and X Signified the displacement of the power piston and the input 
piston valve rod assembly. Ee was the load pressure at the power 


piston, and K. was the fluid spring force at the power piston. 


te 
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Scheduling of Control Parameters. 


Figure 6. 
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B. COMPUTER SIMULATION RESULT 

As mentioned previously, computer simulation results were taken for 
valve stroke displacements, Xa! Peon O,GUS tO O7012 ane Figure 7 
shows the resulting power piston displacement. As can be seen, the sys- 
tem response 1S good. 

Figure 8 portrays the velocity of the power piston. These data 
demonstrate that the velocity of the power eacon TS Ct reculy sp EOror— 
tional to the displacement of the valve input, xX at The natural 
hydraulic frequency of the system was approximately 95 Hz for the con- 
trol unit input servovalve and piston and 110 Hz for the power piston 
assembly. These calculations are shown in Appendix B. 

Figure 9 is a plot of the load pressure, eae for various external 
loads. Without external loading, P was almost zero. After an added 
external load, ee was changed slightly for a range of Fo values from 
Beeo.0 to 40,000.00 lbf. 

Calculations show that the displacement of the power piston upon 
Meaaing 15 0.0016 in. for 40,000.0 lbf. Compared to the ee move- 
ment (10.0 in.), this difference is about 0.016 percent. Thus, the 
displacement of the power piston was almost independent of external 
loading. 

The bulk modulus is a fluid parameter which characterizes the 
compressibility of the liquid. As can be seen from Figure 10, the 
value of the bulk modulus has a slight effect on power piston response 
time. Between bulk modulus values of 100,000 and 300,000 PSI, the 
power piston response was changed from 0.56 to 0.34 sec. Compared to 
a Maximum required time of 19 to 20 sec, during which the power piston 


moved 10.0 in., this difference in power piston response time is very 
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Figure 8. Steady State Velocity Response of Power Piston. 
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Figure 10. Bulk Modulus Effects on Response. 
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small. This is an encouraging result Since minimizing entrapped gas to 
insure high Eo values is a difficult and costly process. 

Figure 11 shows the influence of bulk modulus on power piston 
velocity response. Again, in view of a total system response time 
meaqubprement Of 19 to 20 sec, the initial oscillations that occur within 
the first 0.3 to 0.6 sec of system operation can be considered negligible. 


Finally, in Appendix C 1s shown the computer program developed in 


meas Study and a typical computer output (Figure 15). 
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VI. CONCLUSIONS 


In this study, Several types of C.R.P. systems were analyzed. These 
included: 

1. THE TWO VALVE=PISTON COMBINATION TYPE WITH SOLID CONNECTING ROD 
Results for this system have been presented in the previous 
Sections of this report. 

2. A TWO VALVE-PISTON COMBINATION TYPE WITH CONNECTING FLUID LINE OF 
AC Leculan plee CONSthuction 
b. annular pipe construction 
Ge » double annular pipe construction 

In all of the studies involving fluid lines, the fluid lines had very 
high frequency and very low damping ratio characteristics. Thus, essen- 
tially the influence of the connecting lines on response could be neglec-~ 
ted. Appendix D shows the result of various pipe line response calcula- 
irons. 

Computer simulation proved to be a useful design tool regarding the 
C.R.P. system. Full scale prototype testing, however, is the only con- 
clusive method of proving C.R.P. system response. Unfortunately, testing 
is both costly and time consuming. Thus, the computer simulation 
developed in this study could be used to designate critical testing 
instances and to identify those areas where design efforts would be most 
productive. In this way, the time involved in optimizing system perform- 
ance as well as subSequeént production costs could be reduced. Likewise, 
the present computer simulation could be useful in specifying the correct 


approach for design teSting practices. 
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Pee RNDEX Bs SYSTEM RESPONSE 


The governing equations could be separated into two parts - the input 


part (servovalve andpiston) and the output part (power piston). In this 


Appendix, the hydraulic damping coefficient, f, and the natural frequency, 


W,, Of each part are determined. 


h 


A. INPUT PART (SERVOVALVE AND PISTON) 


From the Figure 13, the characteristic equation is 


where 
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Figure 13. Signal Flow Graph of Input Parts. 
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For nominal values of bulk modulus, the following results 


pertain: 
B, (PSI) w, (HZ) G 
100,000 Soe 17 0.00488 
300,000 159-35 0.00845 
B. OUTPUT PART (POWER PISTON) 


From the Figure 14, the characteristic equation is 


7 2 
WF Xcm, 3? fe PP ) $4 ace on 
M V V 
m g m m 
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Q 
ea 
It 


Mo 


2 2 
(S + R) (S + 20H S +, ) 


Again, for nominal values of bulk modulus, the following results 


were determined: 


B (PST) w, (Hz) i 
100,000 ion 4> 0.00897 
‘ 300,000 191.29 0.04094 





Figure 14. Signal Flow Graph of Output Parts. 
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INITIAL 


CONSTANT 
CONSTANT 


CONSTANT 
CONSTANT 
CONSTANT 
CONSTANT 
CONSTANT 
CONSTANT 
CONSTANT 
CONSTANT 
CONSTANT 
CONSTANT 
CONSTANT 
CONSTANT 
CONSTANT 
DYNAMIC 


AEPENDIX C: 


BE=100000.0 
FF=2000.0 
CEM=1. 54 
QMK=33426. 54 
QAK=2057.08 
XM=69 .8757 
VM=4..9 
VA=749.71 
VC=14326.0 
AM=AN02=0 
AA=57.67 
W=11.938 
PS=400.0 
PSA=600.0 


COMPUTER PROGRAM 


T1=0.0,T2=0.0,T3=0.0,T4=0.0,T5=0.0,T6=0.0 


PARAMETER Z=(0.005,0.01.0.012,0.013) 
ie 

XIN1=Z*STEP(0.0 
XI N2=Z*STEP( 5.0) 
XIN3=Z*STEP( 10.0) 
XING=Z*STEP( 12.5) 
alae ree 
KIN6=Z*STEP( 20.0 
ie leat ea 
XI NS=Z*STEP( 30.0 
XVA=KI N1-XI N2-XIN3+XI NY-XIN5+XI N6+XIN7-XIN8 
YDD=(AM/XM)*PM-( 1./XM)*FC 

MO=UNTCRL(T1, YDD) 
-UNTGHE( bey iD) 
PMD=-(4.*BE*QMK/VC)*Y-(4.*BE*AM/VC)*YD-(4.*BE*CEM/VC)*PM ... 
+(4,*BE*QMK/VC)*X 
PM=I NTGRL(T3,PMD 


FK=0.86*W*( PS-PM)*(X-Y) 


XDD=( FK/VM)*Y-(FK/VM)*X+(AA/VM)*PLA 
XD=INTGRL(T4, XDD) 
K=INTGRL(TS,XD) 


PLAD=-(4.*BE*AA/VA 
PLA=I NTGRL(T6,PLAD 
OUTPUT Y 


PAGE MERGE 
METHOD RKSFX 
TIMER FINTIM=35.0,OUTDEL=0. 5,PRDEL=0. 5, DELT=0.002 


& ND 
STOP 
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PEERNUIx De “FLUID LINE RESPONSE 


When a fluid is used for transmitting pressure variations, it is 
necessary to determine the response of the transmission system. In 
the course of this study, the following line response determinations 


were formulated. 
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Figure 16. Circular Type of Line Response Determinations. 
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Figure 17. Annular Type of Line Response Determinations. 
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Figure 18. Double Annular Type of Line Response Determinations. 
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